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The elliptic azimuthal anisotropy coefficient (v2) of the identified particles at midrapidity (|n| < 0.8) was 


investigated in p—Pb collisions at 


SNN = 5.02 TeV using a multi-phase transport model (AMPT). The calcu- 


lations of differential v2 based on the advanced flow extraction method of light flavor hadrons (pions, kaons, 
protons, and A) in small collision systems were extended to a wider transverse momentum (pr) range of up to 
8 GeV/c for the first time. The string- melting version of the AMPT model provides a good description of the 
measured pr-differential v2 of the mesons but exhibits a slight deviation from the baryon v2. In addition, we ob- 
served the features of mass ordering at low pr and the approximate number of constituent quarks (NCQ) scaled 
at intermediate pr. Moreover, we demonstrate that hadronic rescattering does not have a significant impact 
on v2 in p—Pb collisions for different centrality selections, whereas partonic scattering dominates in generating 
the elliptic anisotropy of the final particles. This study provides further insight into the origin of collective-like 
behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy. 
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I. INTRODUCTION 


The main goal of heavy-ion collisions at ultrarelativistic 
energies is to explore the deconfined state of strongly interact- 
ing matter created at a high energy density and temperature, 
known as the gluon plasma (QGP) [1, 2]. An important ob- 
servation for investigating the transport properties of the QGP 
is anisotropic flow [3, 4], which is quantified by the flow har- 
monic coefficients v, obtained from the Fourier expansion of 
the azimuthal distribution of the produced particles [5, 6]: 
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where ọ is the azimuthal angle of the final-state particle angle 
and W,, is the symmetry plane angle in the collision for the n- 
th harmonic [7, 8]. The second-order coefficient v2, referred 
to as the elliptic flow, is derived from the initial state spa- 
tial anisotropy of the almond-shaped collision overlap region 
that is propagated to the final state momentum space. The 
magnitude of the elliptic flow is sensitive to the fundamental 
transport properties of the fireball, such as the temperature- 
dependent equation of state and the ratio of shear viscosity to 
entropy density (7/s) [9, 10]. 
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Over the past few decades, various measurements of el- 
liptic flow in heavy-ion collisions performed at the relativis- 
tic heavy-ion collider (RHIC) [11-14] and the Large Hadron 
Collider (LHC) [15-18] have helped build a full paradigm of 
the strongly coupled QGP. 


Comprehensive measurements of py-differential elliptic 
flow of the identified particles were conducted by the ALICE 
Collaboration [19, 20]. The observed mass-ordering effect 
(i.e., heavier particles have a smaller elliptic flow than lighter 
particles at the same pr) at low pr is well described by hy- 
drodynamic calculations and is attributed to the radial expan- 
sion of the QGP [21]. At intermediate pr, the grouping of v2 
of mesons and baryons was observed, with mesons exhibiting 
less vg than baryons. These behaviors can be explained by the 
hypothesis that baryons and mesons have different production 
mechanisms through quark coalescence, which has been fur- 
ther investigated using constituent quark (NCQ) scaling [22- 
25]. Interestingly, such flow-like phenomena have been ob- 
served in small-collision systems. Long-range double-ridge 
structures were first measured in high-multiplicity pp and p- 
Pb collisions by the ALICE, ATLAS, and CMS collabora- 
tions [26-28]. The measurement of elliptic and triangular 
azimuthal anisotropies in central 3Het+Au, d+Au, and ptAu 
collisions performed by the STAR Collaboration [29] sug- 
gests that sub-nucleon fluctuations also play an important 
role in influencing the flow coefficients in these small col- 
lision systems. In addition, these measurements were ex- 
tended to the identified particles associated with the discovery 
of a significant positive v2 [30, 31]. The observed particle- 
mass dependence of v2 is similar to that measured in heavy- 
ion collisions [30]; however, the origin of such collective- 
like behavior remains unclear. Several theoretical explana- 
tions relying on either the initial state or final state effects 
have been proposed to understand the origin of azimuthal 
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anisotropies in small systems. Studies that extend hydrody- 
namics from large to small systems based on final-state ef- 
fects can well describe v2 of soft hadrons [32—36]; however, 
they are based on the strong assumption that there is suffi- 
cient scattering among constituents in small systems. Hydro- 
dynamics combined with the linearized Boltzmann transport 
(LBT) model can also describe the identified particle və in 
a high-multiplicity small-collision system at an intermediate 
pr [37]. Color-Glass Condensate (CGC)-based models and 
IP-Glasma models that consider the effect of momentum cor- 
relations in the initial state can quantitatively describe some 
features of collectivity in p—Pb collisions [38, 39], but with- 
out clear conclusions, particularly regarding the dependence 
on collision systems and rapidity. 

In addition, an approach called parton escape shows 
that few scatterings can also create sufficient azimuthal 
anisotropies, which have been investigated using multiphase 
transport (AMPT) [40, 41]. The v2 values of light hadrons 
measured in p—Pb collisions are well described in AMPT, 
where the contribution of anisotropic parton escape rather 
than hydrodynamics plays an important role [41]. In this 
study, we extend the AMPT calculations of the pyp-differential 
Vo for identified particles (n+, K*, p(p), A(A)) to higher pr 
region in p—Pb collisions at 5.02 TeV, in order to systemati- 
cally test the mass-ordering effect and baryon-meson group- 
ing at low- and intermediate-py, respectively. We also inves- 
tigate how the key mechanisms implemented in AMPT, such 
as the parton cascade and hadronic rescattering, affect ellip- 
tic anisotropy in small collision systems. In addition, vari- 
ous nonflow subtraction methods with different sensitivities 
to jet-like correlations were studied. 


Il. MODEL AND METHODOLOGY 


A multiphase transport model 


The string-melting version of the AMPT model (v2.26t9b, 
available online) [40, 42] was employed in this study to cal- 
culate vz of the final-state particles in high-multiplicity p- 
Pb at 5.02 TeV. The AMPT model includes four main pro- 
cesses: initial conditions, partial scattering, hadronisation, 
and hadronic interactions. The initial conditions are gener- 
ated from the heavy ion jet interaction generator (HIJING) 
model [43, 44], where minijet partons and soft-excited strings 
are produced and then converted to primordial hadrons based 
on Lund fragmentation. Under the string-melting mechanism, 
primordial hadrons are converted into partons, a process de- 
termined by their flavor and spin structures. Elastic scatter- 
ing between the partons was simulated using Zhang’s parton 
cascade (ZPC) model [45], which includes two-body scatter- 
ing with a cross-section described by the following simplified 
equation: 
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In this study, the strong coupling constant a, was set to 0.33, 
and the Debye screening mass u = 2.2814 fm", resulting in 
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a total parton scattering cross section of o = 3 mb. To isolate 
the effect of partonic scattering, ø is adjusted to be close to 0 
by increasing u (see set "W/o parton scat." in Tab. 1). Once 
the partonic interaction ceases, hadronization with a quark co- 
alescence model is implemented to combine the nearest two 
(or three) quarks into mesons (or baryons) [40]. The formed 
hadrons enter the subsequent hadronic rescattering process 
using a relativistic transport (ART) model [46], in which 
both elastic and inelastic scattering are considered for baryon- 
baryon, baryon-meson, and meson-meson interactions. The 
hadronic interaction time was set by default to tmar = 30 
fm/c. Alternatively, tmaz is set to 0.4 fm/c to effectively turn 
off the hadron scattering process while still considering the 
resonance decay [47](see set "w/o hadron scat." in Tab. 1). 
In addition, the random orientation of the reaction plane was 
turned on and the shadowing effect was considered in this 
analysis. 


Table 1. Details of three configurations 


Description | o(mb) tmax(fm/c) 
w/ all 3 30 
w/o parton scat. ~0 30 
w/o hadron scat. 3 0.4 


Two-particle correlation and nonflow subtraction 


The two-particle correlation (2PC) method is widely used 
to extract the flow signal in small collision systems because 
it can suppress the non-flow contribution from long-range jet 
correlations [26—28, 30, 48]. Similar to Eq. 1, the azimuthal 
correlation between two emission particles can be represented 
by NPS pairs of emitted particles (labeled as C(Ay)) as a 
function of the relative angle Ay = y*—y? between particles 
a and b and expanded in the Fourier series as follows: 


pair 


x 1+25 Vaa (p$, phcos{n(Ay)], 


n=1 
(3) 
where Vna refers to the two-particle n-th order harmonic. In 
a pure hydrodynamic scenario, because particle emissions are 
independent, V,,(p%, pè) can be factorized into the product 
of a single-particle flow v2 and vt: 


d 
C(Ag) = ake 


Vna (PT, pr) = Vn (PT) Vn (pr). (4) 


Based on the factorization assumption, v,, of a single parti- 
cle a can be obtained using the 3x2PC method, which was 
recently proposed by the PHENIX Collaboration [49]. This 
requires the formation of two-particle correlations between 
three groups of particles (labeled a, b and c) and the extrac- 
tion of the flow coefficients for three combinations: 


N Vna (PS, ph) Vna (PS, p$) 
nla) = y a DPT., 


(5) 
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In small-collision systems, two main types of nonflow con- 
tributions to the flow signal are the near-side jet and away- 
side jet (recoil jet) correlations. The former can be effec- 
tively removed by introducing a large rapidity gap between 
the trigger and associated particles during the construction 
of the correlations. Several methods have been developed 
to subtract the latter [50]. A traditional approach is to di- 
rectly subtract the correlation function distribution obtained 
from low-multiplicity events [27, 30] from that obtained from 
high-multiplicity events. This method assumes that the yield 
and shape of dijets are identical for both collision types as 
follows: 


OM (Aga C™(Ap)x 1+2 5 Vracos|n(Ay)] 


n=1 


=a +2 ` a,cos[n(Ay)], 


n=1 


(6) 


where C!M(Ay) and C™M (Ag) represent the correlation 
function distributions obtained for low- and high-multiplicity 
events, respectively. This method relies on the "zero yield 
at minimum" (ZYAM) hypothesis [27, 30] that a flat combi- 
natoric component should be subtracted from the correlation 
function in low-multiplicity events. Therefore, the fit parame- 
ter a is the absolute modulation in the subtracted correlation 
function distribution and characterizes the modulation rela- 
tive to a baseline, assuming that such a modulation is not 
present in the low-multiplicity class below the baseline. In 
this case, the flow coefficient Vp a is calculated as 


Vna = an/(ao + b), (7) 
where b is the baseline, estimated using the minimum cor- 
relation function for low-multiplicity events. However, the 
measurement of jet-like correlations in p—Pb collisions indi- 
cates that the dependence of the dijet yield on the particle 
multiplicity cannot be ignored. In this case, a new template 
fit method was developed by the ATLAS collaboration [51], 
where the correlation function distribution obtained in high- 
multiplicity events is assumed to result from the superposition 
of the distribution obtained in low-multiplicity events scaled 
up by a multiplicative factor F and a constant modulated by 
cos(nAyw) for n > 1, as shown in 


3 
C(Ag) = FC™ (Ay) +G +2 X Vaacos(nAy)), 


n=1 
(8) 
where G denotes the normalization factor that maintains the 
integral of C(Ay) equal to CHM (Ay). Furthermore, an im- 
proved template fitting method [52] developed in recent years 


was tested. It applies a correction procedure to the default ? 


template fit method by considering the multiplicity depen- 
dence of the remaining ridge in low-multiplicity events, as 
shown in 


FGM 2 2 
Vna = Vna (tmp) zZ gml nA (tmp) E Vaa (LM)), 
(9) 
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where V,,a (tmp) and V,”, (LM) are obtained by using the de- 
fault template method for high- and low-multiplicity events. 
All these nonflow subtraction methods are implemented in 
this study, and their different sensitivities to nonflow effect 
are also discussed. 


HI. ANALYSIS PROCEDURES 


To directly compare the AMPT calculations with the re- 
sults from ALICE, we focused on the particles within the 
pseudorapidity range |7| < 0.8, aligning with the TPC ac- 
ceptance in ALICE [53]. In the 3x2PC method, long-range 
correlations were constructed between the charged particles 
at mid-rapidity, forward rapidity (2.5 < ņ < 4), and back- 
ward rapidity (—4 < 7 < —2.5), that is, the central-forward 
correlation (—4.8 < An < —1.7), central-backward correla- 
tion (1.7 < An < 4.8), and backward-forward correlations 
(—8 < An < —5). In addition, the centrality classes are de- 
fined by counting the charged particles in the acceptance of 
the VOA detector [53], that is, 2.8 < 7 < 5.1. 

The correlation function distribution C(Ay) was obtained 
by correcting the number of particle pairs in the same events 
normalized to the number of trigger particles Ntrig by using 
an event-mixing technique: 

1 d? Npairs S(Ag, An) 
Nirig dAndApy  B(Ay, An)’ 


C(Ag, An) = (10) 


where S(Ay, An) = ++ d Nsame ig the correlation func- 


Nurig dAndAg 

tion in same events and B(Ay, An) = ose is the as- 
sociated yield as a function of Ay and Ay in mixed events. 
Factor a is used to normalize B(Ay, An) to unity in the Ay 
region of the maximal pair acceptance. The obtained 2-D cor- 
relation function C(Ay, An) is projected onto Ay axis, and 
we follow the nonflow subtraction procedures and factoriza- 
tion, as discussed in Eq. 5—Eq. 9, v2 of the charged particles 
at |7| < 0.8 can be calculated. 


IV. RESULTS AND DISCUSSIONS 


We first investigated the pr spectrum of the identified par- 
ticles before performing the flow analysis. Figure | illustrates 
the pr distribution of proton, pion, and kaon in 0-20% high- 
multiplicity p—Pb collisions at \/snn = 5.02 TeV, which are 
obtained from AMPT with three different sets of configura- 
tions listed in Tab. 1 and ALICE experimental data [54]. 
The AMPT results, both with and without hadronic rescat- 
tering, are consistent. This behavior differs from previous 
findings in heavy-ion collisions, where the hadronic interac- 
tion significantly reduces the particle yield [55]. The spec- 
trum obtained in the AMPT without considering the parton 
cascade process is enhanced compared to that obtained with 
partonic scattering, and this enhancement is more significant 
at a high pr. This outcome is expected because partons expe- 
rience energy loss during the parton cascade, which reduces 
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Fig. 1. (Color online) The pr distribution of pions, kaons, and pro- 
tons in 0-20% high-multiplicity p—Pb collisions at \/snn = 5.02 
TeV, obtained from AMPT model calculations, is compared to AL- 
ICE measurement [54]. The results in AMPT without hadronic scat- 
tering and partonic scattering are also presented. 


the production of final-state particles. In addition, the ratios 
of the pr spectra obtained from the AMPT calculations and 
data are shown. The AMPT model calculation reproduces the 
particle yields well at low and intermediate py values when 
both partonic and hadronic scattering are included; however, 
it overestimates the high pr data because parton-parton in- 
elastic collisions and, subsequently, hard parton fragmenta- 
tion are absent in the model [42]. 


Figure 2 (left) shows the v2 of pions, kaons, protons and 
A as a function of pr in 0-20% high-multiplicity p—Pb col- 
lisions at \/snn = 5.02 TeV, obtained in AMPT calculations 
with 3x2PC method. A comparison with the ALICE mea- 
surement for vg of charged hadrons, pions, kaons, and pro- 
tons [30] and the CMS measurement for v2 of K? and A [31] 
is also presented. The AMPT calculations applied the tem- 
plate fit method to suppress the away-side jet contribution and 
considered the ZYAM assumption to enable direct compari- 
son with the observed data. The v2 values of charged hadrons, 
pions, and kaons can be described well by AMPT calcula- 
tions, but the v2 values of baryons (protons and A) cannot be 
reproduced. In addition, the mass-ordering effect (i.e., the 
vg of baryons is lower than that of mesons) is reproduced 
for pr < 2 GeV/c. Owing to the advanced flow extraction 
method, the calculations of va were extended to the high-pr 
region, up to 8 GeV/c in the AMPT model for the first time. 
The və values of protons and A are consistent, and both of 
them are observed to have a higher value of və value than 
that of the mesons for 2 < pr < 7 GeV/c. The observed 
meson-baryon particle type grouping in heavy-ion collision 
flow measurements indicates collective behavior at the par- 
tonic level, leading to the coalescence of quarks into hadrons. 
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The number of constituent quarks (NCQ) scaling techniques 
described in [22] can be used for further studies of this group- 
ing. v2 and pr in Fig. 2 (left) are replaced by v2/n, and 
pr/Ng, Where the ng is the number of constituent quark in 
mesons (Nng = 2) and baryons (ng = 3), as shown in Fig. 
2 (right). v2/ng obtained from the data show approximate 
values at intermediate py; however, the results calculated in 
AMPT cannot reproduce the scaling in pp/ng > 1 GeV/c. In 
order to consider the observed mass hierarchy of v2, we also 
plot the v2 of identified particle as a function of the transverse 
kinetic energy kEp (kEr = mp — Mo = y på + m2 — mo), 
and its NCQ scaling in Fig. 3 (left), and Fig. 3 (right)). 
All particle species showed a set of similar v values af- 
ter NCQ scaling in kEr/n, < 1 GeV, confirming that the 
quark degree of freedom in flowing matter can also be probed 
in the transport model. However, this NCQ scaling is vio- 
lated for KE; /n, > 1 GeV. This may be attributed to the 
hadronization mechanism implemented in the AMPT model 
used in this study, where baryons are produced only after the 
formation of mesons by simply combining the three nearest 
partons, regardless of the relative momentum among the co- 
alescing partons. This results in an underestimation of the 
baryon və at intermediate pr in this study. An improved co- 
alescence model implemented in the newer AMPT [56] in- 
troduced a new coalescence parameter to control the relative 
probability of a quark forming a baryon instead of a meson 
precisely. This improvement could have different NCQ scal- 
ing on və but requires more systematic studies. Further stud- 
ies on v2 calculations in small collision systems with other 
improved hadronization mechanisms, for example, consid- 
ering the Wigner function [57] and hard parton fragmenta- 
tion [58], should be performed in the future. 


We also extend our investigation to include a study of inte- 
grated v2 within various centrality bins spanning the 0-60% 
range. We focus on the region where the NCQ scaling cri- 
terion is satisfied, that is, for transverse kinetic energies per 
constituent quark (kE7;/n,) ranging from 0.4 to 1 GeV. The 
non-flow contribution was estimated and subtracted within 
the 60-100% centrality class by using the template fit method. 
As shown in Fig. 4, the v2 values as a function of centrality 
exhibit a systematic decrease from central to peripheral colli- 
sions, reflecting the changing dynamic conditions and particle 
production mechanisms in different collision zones. Intrigu- 
ingly, in the v2 measurements, we observed a distinct mass- 
splitting phenomenon, with baryons and mesons exhibiting 
distinct elliptic flow patterns. Such a mass dependence in v2 
is similar to that in heavy-ion collisions at the LHC energies 
presented in a previous study [47]. This provides valuable in- 
sights into the collective behavior of different particle species 
within the evolving fireball created during these collisions. 


Moreover, to gain a deeper understanding of the NCQ scal- 
ing properties, we explored the ratios of ng-scaled integrated 
və values for protons relative to pions and kaons relative to 
pions as functions of centrality. The results are shown in Fig. 
5. A notable trend is observed in these ratios: they tend to ap- 
proach unity as the collisions become more peripheral. It in- 
dicates that the collective flow of particles in low-multiplicity 
events may be approaching a behavior that is closer to the ex- 
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Fig. 2. (Color online) Left: the v2 as a function of pr in 0-20% high-multiplicity p—Pb collisions at ,/snn = 5.02 TeV, obtained from default 
AMPT model calculations with 3x 2PC method, is compared to ALICE and CMS measurement [30, 31]. Right: the pr-differential v2 scaled 


by the number of constituent quark (nq). 
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Fig. 3. (Color online) Left: the v2 as a function of transverse kinetic energy (KEr) in 0-20% high-multiplicity p—Pb collisions at ,/snn = 
5.02 TeV, obtained from default AMPT model calculations with 3x2PC method, is compared to ALICE and CMS measurement [30, 31]. 
Right: the kEr-differential v2 scaled by the number of constituent quark (n4). 
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Fig. 4. (Color online) The integrated v2 in 0.4< kEr/ng <1 GeV 
for pion, kaon and proton varying with the centrality. 
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The effects of partonic and hadronic scattering on the el- 
liptical anisotropy of the final-state particles were examined 
in this study. Figure 6 shows the calculated py-differential 
v2 of pions, kaons, and protons in AMPT with and without 
considering hadronic rescattering process in 0—20% high- 
multiplicity p—Pb collisions. The results show that the ra- 
tio of the vo values with and without hadronic rescattering is 
consistent with unity for all particle species, indicating that 
the hadronic rescattering mechanism has almost no effect on 
v2 in high-multiplicity p—Pb collisions. We also investigated 
the centrality dependence of the hadronic rescattering effects 
by calculating py-integrated v2 in several centrality bins be- 
tween 0 and 60%, as illustrated in Fig. 7. The results demon- 
strate that the influence of hadronic rescattering is indepen- 
dent of the centrality selection and has almost no impact on 
NCQ scaling in the range of 0.4< kEr/ng <1 GeV. 

On the other hand, when we set the parton scattering cross- 
section o to zero but maintain the hadronic scatterings, the 
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Fig. 5. (Color online) The ratio of integrated v2 within 0.4< 
kEr/nq <1 GeV for proton over pion and kaon over pion vary- 
ing with the centrality. The dash line represents the location of unity 
ratio. 
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Fig. 6. (Color online) The pr-differential v2 of pions, kaons, and 
protons calculated in AMPT model with and without considering 
hadronic scattering. The ratios of the two sets are also presented. 


central-backward (CB) correlations is almost zero, as shown 
in Fig. 8. If both the partonic and hadronic scatterings are 
turned off, the results remain consistent with zero. This indi- 
cates that the elliptical anisotropy in high-multiplicity small- 
collision systems is mostly generated by parton scattering. 
Our conclusion is consistent with previous studies on the 
AMPT [41], which suggested that the majority of elliptic 
anisotropies comes from the anisotropic escape probability 
of partons. 

Finally, different non-flow subtraction methods were in- 
vestigated in this study. Figure 9 (left) shows the pr- 
differential v2 of the charged particles calculated using the 
3x2PC method in 0—20% high-multiplicity p—Pb colli- 
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Fig. 7. (Color online) The integrated v2 in 0.4< kEr/ng <1 GeV 
for pions, kaons, and protons calculated in AMPT model with and 
without considering hadronic scattering. The ratios of the two sets 
are also presented. 
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Fig. 8. (Color online) The pr-differential V2, for central-forward 
(CF) and central-backward (CB) correlations calculated in AMPT 
model with and without considering partonic scattering. 


sions. Several nonflow subtraction methods are implemented. 
To demonstrate how the nonflow contribution is removed, v2 
obtained with a direct Fourier transform of the C(Ay) cor- 
relation (as shown in Eq. 3). The results show significant 
suppression across all the subtraction methods, particularly at 
higher pr values where jet correlations are dominant. The re- 
sults obtained with peripheral subtraction and template fitting 
were consistent, indicating that the away-side jet contribution 
was automatically removed using the 3x2PC method, even 
though the dependence of the jet correlation on multiplic- 
ity was not considered in the peripheral subtraction method. 
The və calculated using the improved template fit method was 
slightly lower than that from the template fit, and it was simi- 
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Fig. 9. (Color online) The pr-differential v2 of charged hadrons 
calculated in AMPT with different nonflow subtraction methods. 


377 lar to the features observed in the ATLAS measurement [52]. 
373 The same conclusions were drawn for the extraction of the 
379 identified particles (pions, kaons, protons, and A) və. 


V. SUMMARY 
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This study systematically investigated the elliptic 
382 anisotropy of identified particles (pions, kaons, protons, and 
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